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Abstract 

We replace time-averaged entanglement by ensemble-averaged entan- 
glement and derive a simple expression for the latter. We show how to 
calculate the ensemble average for a two-spin system and for the Jaynes- 
Cummings model. In both cases the time-dependent entanglement is 
known as well so that one can verify that the time average coincides with 
the ensemble average. 

1 Introduction 

The entanglement of particles is in principle a time-dependent quantity. This 
time-dependence has been analysed recently in chaotic systems [1 , in exper- 
imental spectra of triatomic molecules [2J, and in Rydberg atoms [3J. Time- 
dependent entanglement has been studied in theoretical models, like the Dicke 
model [I], a model of coupled kicked tops [5], the Harper Hamiltonian [B], a 
dimer model [7], Bose-Einstein condensates [§]. In these papers, the notions 
of time-averaged entanglement and ensemble-averaged entanglement have been 
shown to be useful in monitoring phase transitions, although the generality of 
this relationship has been questioned, see e.g. [3j. 

In addition, time-averaged and ensemble-averaged entanglement are con- 
served quantities of quantum microcanonical ensembles [9j. As such they are 
of interest in the study of closed systems. This context is suited to discuss the 
relation between both concepts, and is the starting point of the present paper. 

The entanglement of formation of a pure state is taken here to be defined 
as the von Neumann entropy of the reduced density matrix. Often, the von 
Neumann entropy is replaced by the linear entropy because the latter can be 
computed more easily. In the present paper the use of the linear entropy is 
essential to obtain simple results. 

For the sake of completeness, and to fix notations, the definitions of entan- 
glement of pure and of mixed states are reproduced in the next section. Section 
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3 introduces the entanglement of microcanonical ensembles of wavefunctions. 
The main result is announced in Section 4. The proof is found in the Appendix. 
It is followed by a section devoted to the complications that arise when the state 
of the system has additional symmetries. In Sections 6 and 7 the main result is 
applied to a system of two interacting spins. Section 8 deals with the Jaynes- 
Cummings model. For this model the time-dependent entanglement is known so 
that it can be compared with the ensemble average. A short discussion follows 
in Section 9. 



2 Definition of entanglement 

Consider two independent subsystems, labelled A and B. With each normalised 
wavefunction ip of the combined system corresponds a reduced density matrix 
p A of the subsystem labelled A. The latter is defined by the relation 

Tr A p A X = (ip\X <g) I\ip) for all X. (1) 

The entanglement of ip is then equal to the von Neumann entropy of pa- For 
technical reasons we replace this entropy by the linear entropy. The general 
definition of entanglement is 

£ A (iP) = S f (p A )=Tr A p A f(p A ), (2) 

with f(x) = — lnx in the von Neumann case, and f(x) — 1 — x in the case of 
the linear entropy. 

If ip is of the product form ip = ip A ® ip B then p A is a one-dimensional 
projection operator. Hence, the entanglement vanishes. A similar definition 
holds for p B and for £ B . The entanglements £ A (ip) and £ B (ip) are equal [ID] . 
To see this, select a basis u m in subsystem A and a basis v p in subsystem B, so 
that 

= X! VPn U n ® V n . (3) 
n 

This is possible by means of the Schmidt construction. Then p A and p B are 
diagonal, with eigenvalues p n , and with entropy equal to ^2 n Pnf{pn)- 

Often, the state of the system is not described by a wavefunction ip but by 
a density matrix p. Such a density matrix can be written into the form 

P = ^2Pn\lpn) (tpn\, (4) 
n 

with p n > 0, ^2 n Pn — 1, and with ip n normalised wavefunctions. Then the 
entanglement of p has been defined [11] as the minimum of the average entan- 
glement 

£(p) = min^2p n £ A (ip n ), (5) 

n 

where the minimum is taken over all possible ways to write ([4]). 
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3 Definition of mean entanglement of a micro- 
canonical ensemble 



The mean entanglement, which is studied in the present paper, is not the average 
l[5]). but rather the average over a microcanonical ensemble, as introduced in [5]. 

Let be given a density matrix p, which is diagonal in the orthonormal basis 
of wavefunctions ip n , with eigenvalues p n : pip n = p n ip n - Associated with this 
diagonal density matrix is an ensemble of wavefunctions of the form 

tp = ^2VKe iXn ip n , (6) 

n 

where the Xn are arbitrary phase factors. The ensemble average of the entan- 
glement is then denoted £ and is given by 



£ = EvS"H ) , (7) 




where the average over \ is obtained by integrating over all phase factors \ n 
from to 2ir, normalised by dividing by a factor 2ir. Note that the ensemble 
average ([7} does not depend on the chosen subsystem because £ A {fa) = £ B (i[>) 
for all ip. 

The ensemble ([6]) can be obtained by starting from a single wavefunction ip, 
in combination with the quantum mechanical time evolution. The Hamiltonian 
H is the generator of the unitary time evolution 

fa = U(t)ip with U{t) = e- lh ~ ltH . (8) 
The time average of the entanglement £ A (ip) is then defined by 

(£A4>t)) t = lim i f T dt£ A (fa). (9) 

Assume now that the Hamiltonian is diagonal in the basis of wavefunctions ip n , 
with eigenvalues e n . Then one has 

^ t = J2 x " e ~ m ~ lent ^ ( 10 ) 

n 

and hence 

{£ A ^t)) t =Km o ±£dt£ A ^A„e-^ le "*^ ■ (H) 

The wavefunctions fa belong to the ensemble ([6]) with p n — |A„| 2 . If the condi- 
tions of the classical ergodic theorem holds then the time average (fTTj) coincides 
with the ensemble average (0 — see [5]. But even when the classical ergodic 
theorem does not hold one can continue to use the ensemble average instead of 
the time average because experimentally the slightest perturbation may restore 
ergodicity. 
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4 Main result 



Let be given an ensemble of wavefunctions of the form ((6|) . With each of the 
basis vectors ip n is associated a couple of reduced density matrices p A and p B . 
For convenience, these will be denoted a n and r„. Introduce the density matrices 
a and t, defined by 

<T = ^p„er„ and r = ^p„r„. (12) 

n n 

In some sense, a is the ensemble average of p A1 r is the ensemble average of p B . 

Our main result is now that the mean entanglement, defined by ([7]), using 
the linear entropy Si, can be written as 

£ = Si(a) + Si(t) - A, (13) 
where A is a contribution common to both Si(a) and Si(t). It is given by 

A = 1 - Pi Tr A o* = 1 - £ Tr B r, 2 „. (14) 

m m 

The proof of this relation is given in the Appendix. 

The applications of (fTB")) are explored in later sections. 

5 Degeneracies 

The ensemble (j6|) is uniquely defined by the density operator p in the case that 
the eigenvalues p n of p are two-by-two distinct. Then the eigenfunctions ij) n are 
unique up to a phase factor. However, If some of the eigenvalues p n coincide 
then the the orthonormal basis is non-unique. In particular, if a non-zero eigen- 
value p n is degenerate then different choices of orthonormal wavefunctions may 
influence the value of £. This shows that £ is the average entanglement of the 
ensemble and is not suitable as a definition of the entanglement of p. 

A similar question is whether the entanglement of the ensemble can be used 
as the definition of the mean entanglement of the wavefunction ip. Consider 
the situation that some of the eigenvalues e„ of the Hamiltonian H are de- 
generate. Then the basis of wavefunctions tf> n , which diagonalises H, is not 
uniquely defined (up to phase factors). In that case the wavefunction ip should 
not be decomposed into an arbitrary diagonalising orthonormal basis. Rather, 
it should be projected onto the invariant subspaces of H. This determines in a 
unique way an orthonormal basis which then can be used to form the ensemble 
associated with ip. An example of the degenerate case follows below. 
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6 Two-spin example 



The simplest example is that of two quantum spins, each described by Pauli spin 
matrices, and a Hamiltonian H which is diagonal in the basis of wavefunctions 

= ^(ITT> + IU» 

4>2 = -±=(|TT)-|U» 

^3 = ^(ITD + IIT)) 

^4 = ^(lU)-IIT)). (15) 

We assume that the energy levels are non-degenerate. Their actual value is not 
needed. 

The reduced density matrices a m and r m , corresponding with tjj mi are all 
equal to il. Hence, also the averages a and r are equal to |l. As a consequence, 
the linear entropies Si(a) and Si (r) both equal 1/2. However, the common part 
A depends on the choice of wavefunction ip. One finds 

A = 1 - E Prn Tr A =l~\Y,tin- ( 16 ) 

m m 

The final result for the mean entanglement of is therefore 

\ m / m 

Note that this result lies between 1/4 and 1/2. 

It is possible but tedious to verify by explicit calculation that the mean 
entanglement (fl7|) coincides with the time average of £ A (ipt), as it should be. 



7 Degenerate two-spin example 

Consider a two-spin system with energy — e for anti-parallel spins and +e for 
parallel spins. This is a degenerate limit of the previous example. The ensemble, 
generated by an arbitrary ip, now contains two free phase factors instead of four. 
It consists of all wavefunctions of the form 

e i * 1 \ 1 ip 1 + e ix *\ 2 ip2 (18) 
with xi and \2 arbitrary, and with \\ipi, ^2^2 given by 

= I Tl) (Tl \i>) + 1 IT) (IT \i>) 
A 2 V 2 = I TT) (TT W) + 1 II) (II IV>>. (19) 
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The coefficients Ai and A2 are chosen in such a way that ip% and ip2 are nor- 
malised. The reduced density matrices are found to be 



P1&1 = 



P+- 
p_+ 



>** - ( T ) ■ < 20 > 

with p_| = KtJ. \ip) \ 2 and similar definitions for p |_, p ++ and p Similar 

expressions hold for t\ and ti 



P\T\ 



P-+ 
p+_ 



M = ( n ++ „ ° ) • (21) 



p 

The mean entanglement can now be calculated using (|13p . The result is 

£ = 2 P++ p__+2 P+ _p_ + . (22) 
In the notation of the previous section this becomes 

£ = \{Pi-P2? + \{P3-Pa) 2 , (23) 
which is less than (|17p by the term —p\P2 — PsPa- 



8 The Jaynes-Cummings model 

The Jaynes-Cummings model [121 113] describes a two-level system interacting 
with a harmonic oscillator. The latter represents a single mode of the electro- 
magnetic field in a cavity. The model has been studied extensively. 
The Hamiltonian of the model reads 

H = huja) a + — hhjQO z + hK(a) O- + aa + ), (24) 

with and a creation and annihilation operators of the harmonic oscillator, 
and with the Pauli matrices a z and a± describing the two-level system. 

Let \g) and |e) denote the ground state, respectively the excited state of the 
two-level system. Let |n), n = 0, 1, • • • denote the eigenstates of the harmonic 
oscillator. The eigenstates of the Jaynes-Cummings Hamiltonian are explicitly 
known, see e.g. |14j . An orthonormal basis of eigenfunctions is given by 

V-o = |fl)<8|0) 
V>i,„ = cos(0 n )|ff) ® + +sin(<9 n )|e) ® \n) 
i>2,n = -sm(9 n )\g) ® |n+ 1) +cos(6»„)|e) ® \n). (25) 



G 



The angles 6 n , n = 0, 1, 2, ■ • • follow from the relation 



tan0 n = Ky- U+ 1 (26) 

5 (w - u> ) + A„ 



with 



A„ = yi(^-^ ) 2 + K 2 (n+l). (27) 

The time-dependence of the reduced density matrix can be calculated ex- 
plicitly if the initial state is a product state with the two-level system in the 
excited state and the harmonic oscillator is in the n-th eigenstate, see e.g. [T3]. 
The result for the reduced state of the two-level system is 

p A (t) = W n (t)\g) (g\ + (1 - W n (t))\e) (e\ (28) 

with 



W n [t) = 2 7 „ sin 2 (tKy/n + l) (29) 

and with 

7n = isin 2 (20„). (30) 

The linear entanglement is therefore 

£ A ^ t ) = l-Wl{t)-{l-W n {t)) 2 . (31) 
The time average equals 

(£ A (^)) t = 2 7n -3<£ (32) 

It is straightforward to calculate the reduced density matrices for the eigen- 
functions of the model. The result is 

(33) 
(34) 
(35) 
(36) 
(37) 
(38) 

Hence, it is straightforward to evaluate the mean entanglement £ for an arbitrary 
wavefunction tp. However, we did not succeed to rewrite the resulting expression 
in a simple and transparent way. 

In the case that ip is of the product form 

%j) = |e) ® |n) = sin(6' ri )'0i i „ + cos(6> n )V>2,n (39) 



00 


= WW, 






= cos 2 6»„n+l) 


(n + l\ + sin 2 9 n n) (n 


0"2,n 


= sin 2 9 n \n + 1) 


(n + l\+ cos 2 9 n \n) (n 


TO 


= \g)(g\, 




Tl, n 


= cos 2 9 n \g) (g\ - 


fsin 2 9 n \e) (e|, 


T2,n 


= sin 2 6»„ \g) (g\ - 


1- cos 2 6»„ e) (e . 
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one obtains 



a 



sin 2 (0 n )<7i,„ + cos 2 (6>„)er 2 ,ri 
7„|n+l)(n + l| + (l-7„)|n) (n| 



T 



sin 2 (6' n )ri,„ + cos 2 (6> n )r 2 , 
ln\9) (g\ + (1 - 7n) |e) (e| 



A 



1 - sin 4 (#„) Tr B r 2 „ - cos 4 (0„) Tr B r| 
l-(l-7«) 2 . 



(40) 



This leads to the result 



£ 



27n - H 



2 



(41) 



which is identical with the time-averaged result (132)) . 

9 Discussion 

The calculation of the time-dependence of the entanglement of a quantum sys- 
tem is a hard problem. The average over time is more accessible because it can 
be replaced by an ensemble average. This is in particular so when the entangle- 
ment is defined using the linear entropy instead of the von Neumann entropy 
of the reduced density matrix, because in that case there exists a simple ex- 
pression for the mean entanglement — see (1131) . We have used this expression 
in a two-spin system and in the Jaynes-Cummings model. For these systems 
it is feasible to calculate both the averages over time and over the ensemble of 
wavefunctions. The results of the two calculations coincide, as it should be. 

We have pointed out that a systematic degeneracy of the energy levels of 
the Hamiltonian due to the presence of a symmetry influences the choice of the 
ensemble of wavefunctions, used in the calculation of the average entanglement. 
In the example of the two spin system the additional symmetry leads to a 
reduction of the entanglement. 

Finally let us note that the use of the linear entropy is rather essential in 
our paper. It is of course possible as well to define the ensemble average of the 
entanglement, based on the von Neumann entropy. It is however unlikely that 
a simple formula like (fl3)) exists in that case. 

Appendix 

Here, the proof of (|T5|) is given. From the definition ([7]) follows, assuming a 
linear entropy, 



£ 



l-(Tr A p 2 ( X )) x , 



(42) 



with p A (x) defined by 
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= ^V&e-^-^'Tr^™! (43) 

mn 

with cr mn defined by 

Tr A a mn X = {^ m \X®I\ilj n ). (44) 

Hence, one obtains 

\ 2v 



p— i{Xm~Xn) n 
). c u mn 



£ = 1 - ( Tr A I ^2 VPmPnt 
\ \ mn 

= 1 - ^ ^PmPnPrPs (e~' (x "' ~ X " } e~' (Xr ~ X ° } \ Tl" A <7 mn (7 rs 
mn rs ^ 

= 1 - E *™P» Ti a <^„„ - £ Tr A a mn a nm 

mn mn 

+ E^ Tr A 17 ™- (45) 

m 

Now use that a mm = a m to see that the first two terms yield the contribution 
Si((t). The last term is 1 — A. Hence, it remains to be shown that 

Si(t) = 1 - £ p m p n Tr A <T mn a nm . (46) 

mn 

Choose an orthonormal basis u r for the subsystem A and an orthonormal 
basis v p for the subsystem B. Then one has 

Umn = £ [(lprn\u s <g> Vp) (u r <g) V p \lj} n )] \r) (s\ (47) 
rsp 

so that 

Tr A a mn a nm = £ {ip m \u s ® v p ) (u r ® u p |^„) 

X(^„|Wr <8) V q ) (u s (g) Wgl^m) 

K) (u g |) (V'nl (I® K) <Up|) IV 1 ™) 



= Tr B T m T n . (48) 
The relation (|4"6"| now follows readily. 



References 

[1] X. Wang, S. Ghose, B.C. Sanders, and B. Hu, Phys. Rev. E 70, 016217 
(2004). 



[2] X.W. Hou, M.F. Wan, and Z.Q. Ma, Chem. Phys. Lett. 426, 469 (2006). 

[3] M. Lombardi and A. Matzkin, Phys. Rev. A 73, 062335 (2006). 

[4] X.W. Hou and B.B. Hu, Phys. Rev A 69, 042110 (2004). 

[5] R. Demkowicz-Dobrzanski and M. Kus, Phys. Rev. E 70, 066216 (2004). 

[6] A. Lakshminarayan and V. Subrahmanyam, Phys. Rev. A 67, 052304 
(2005). 

[7] X.W. Hou, J.H. Chen, and B. Hu, Phys. Rev. A 71, 034302 (2005). 

[8] Q. Xic and W. Hai, Eur. Phys. J. D 33, 265 (2005). 

[9] J. Naudts and E. Van der Straeten, J. Stat. Mech. P06015 (2006), 



arXiv:quant-ph/0602039, 



[10] C.H. Bennett, H.J. Bernstein, S. Popescu, and B. Schumacher, Phys. Rev. 
A 53(4), 2046 (1996). 

[11] C. H. Bennett, D. P. DiVincenzo, J. Smolin, and W. K. Wootters, Phys. 
Rev. A 54, 3824 (1996). 

[12] E.T. Jaynes and F.W. Cummings, Proc. IEEE 51, 89 (1963). 

[13] F.W. Cummings, Phys. Rev. A 140, 1051 (1965). 

[14] A.K. Rajagopal, K.L. Jensen, and F.W. Cummings, Phys. Lett. A 259, 
285290 (1999). 



10 



